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We argue that the classical form of the dipole-dipole interaction energy cannot be used to model 
the interaction of the bosons in a dilute Bose-Einstein condensate made of polar atoms. This fact 
is due to convergence of integrals, if no additional restrictions are introduced. The problem can 
be regularized, in particular, by introducing a hard sphere model. As an example we propose 
a regularization consistent with the long range behavior of the effective potential and with the 
scattering amplitude of the fast particles. 



The first clear observation of pure Bose-Einstein con- 
densation was done using a ultracold gas of neutral 
bosonic atoms [Q. After these remarkable experiments 
an explosive growth of the interest on the subject hap- 
pened [||. Following the experimental guidelines, the 
theoretical analysis of the problem concentrated on con- 
densates made of alkaline atoms (neutral and nonpolar) 
interacting coUisionally through s-wave spherically sym- 
metric effective potentials, which lead to an interaction 
proportional to the local atomic density. 

Although this description of the problem captures most 
of what has been done experimentally up to now, very 
recently there has been increased interest on other, more 
complex, interactions. Specifically, there has been some 
interest on nonlocal interactions (see e.g. and refer- 
ences therein). In fact any realistic interaction should 
always be nonlocal due to the fact that the range of inter- 
action is not zero. Another situation corresponds to non- 
symmetric and nonlocal interactions considered in Refs. 
iQ-Q, such as those which appear when the bosons in- 
side the condensate are polar molecules (with permanent 
or externally induced dipole moments). The possibility 
of finding condensates with nonsymmetric interactions is 
very interesting since it would lead to a bunch of new 
phenomena and many possibilities for control. 

In all the previously cited cases the model for the dy- 
namics of the system is the zero temperature mean field 
theory, i. e. a Gross-Pitaevskii equation, for the conden- 
sate dynamics 
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where 5* is the condensate wavefunction normalized to 
the total number of particles, 
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a is the s-wave scattering length and F is a real function 
taking into account the nonlocal part of the interactions. 
Of course, V vanishes or is very small when the interac- 
tions are purely local. When the interactions are medi- 
ated by dipoles, they are supposed to be long range. In 
such situation it is proposed that the interaction is ruled 
by a potential of the type P,™?! 
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where Gr = (r — r')/|r — r'| and d(r) is the dipole moment 
at point r. The simplest version of Eq. (^), where d(r) = 
d is homogeneous and the potential can be considered as 
a function of relative distances only, has the form 
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where cos 9 = ed • e^. In some cases [Q, instead of the 
simple angular dependence f — 3 cos^ 6, which is propor- 
tional to the spherical harmonic 1^20 more complex depen- 
dencies in the form of a finite combination of spherical 
harmonics may appear. 

It is our intention in this brief report to show that 
all those potentials lead to a singular (in mathematical 
sense) interaction which must be regularized to give phys- 
ically acceptable results. This situation is not surprising 
if one takes into account that even in classical theory (i) 
Eq. (^) is an asymptotics of the dipole-dipole interac- 
tions (see e.g. Q) and thus fails in the near zone, and 
(ii) Eq. (H) displays the realistic dipole-dipole interac- 
tion while V{r) in (|f|) is an effective potential obtained 
within the mean field theory. 

To simplify the analysis, let us first consider the solu- 
tion of the Gross-Pitaevskii equation with an interaction 



f 



term of the form Then the nonhnear nonlocal term 
is proportfonal to the integral 
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being J(r') — |5'(r')|^. Taking into account that due 
to the presence of the trap and the finite norm of the 
solution the wave function decays as |r| — > oo, the only 
point where the integral may lead to a singular behavior 
is r' = r, thus we will employ spherical coordinates and 
estimate the integral in the e— vicinity of u = r — r'. To 
do so we split (hereafter dO. = sin 9d9d(j)) 
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The singular behavior, if present, may only be due to the 
contribution of /e(r) for which we have the following set 
of estimates 
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Taking into account that J(r) does not depend on the 
integration variable one gets an indetermination of type 
X oo which cannot be avoided because the order of in- 
tegration is not defined (and should not be relevant in a 
physically meaningful model) . This is the first indication 
that there may be problems with the type of nonlocal 
interaction kernels considered. 

To rule out the possibility that the singular behavior 
previously found is a product of the specific orientation of 
the dipolcs chosen (i.e. all of them parallel) we consider 
now a general situation. Let us first define a new function 

A(r,u) = d(r + u) -d(r). (8) 

With this definition wc may write 
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Let us formally split this integral into two parts, /(r) — 
/o(r) + /i(r). The first one 
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Then we may write 
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Thus if (|12| ) is convergent one may ensure convergence 
of (^fj) or equivalently (13). Assuming that d(r) is a 
differentiable function and expanding in Taylor series we 
find that A(r,u) = 0{u) and thus the singular behavior 
for small u values is avoided. Thus the convergence of the 
integral depends on the contribution of /o (r) . This means 
that the consideration in what follows can be restricted 
to the model (^ without restriction of generality. 

As we have commented above the integral /o(r) is not 
well defined. Let us show, however that the problem 
may be regularized. Indeed, assuming that the dipole- 
dipole interaction is cut off at some distance Oc one has 
to substitute the lower limit in (^) by Oc what will lead to 
the result of integration which is not divergent at Oc 0. 
Mathematically 

Ioir) = d\r) — dnj{r + u}[l-3cos^9] (14) 

is uniformily bounded with respect to r. The contribu- 
tion of the small scale of the integral can be estimated as 
before by expanding J{r) as 
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In fact one may easily prove that the result may be ex- 
panded in power series of a as 
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so that, in the limit Oc — > one gets a finite value for 
loir). 

The fact that the integral is well behaved when a cutoff 
is used is the reason why no divergences were observed 
in numerical studies |^-^ . The use of any computational 
mesh introduces a cutoff related to the mesh size which 
is effectively equivalent to introducing the hard sphere 
radius. 

The cut-off model, although the most straightforward 
way of avoiding singularities of the problem is not the 
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most natural one, from the physical point of view. In- 
deed, this is not properly a "hard sphere" morel which is 
often used in physics, since the respective cut-off is not 
introduced in the self-interaction term. Thus such a cut- 
off corresponds to zero dipolc-dipolc interaction at small 
distances which is not correct from the physical point of 
view. 

In order to propose another way of regularization in the 
case at hand we take into account that in the field theory 
dipoles appear as an artifact of the energy expansion with 
respect to a small parameter r^jr where r^i is a charac- 
teristic dimension of the dipole (for the next terms of the 
expansion one has to compute multipolar terms). The 
complete potential energy is just made of two-particle 
Coulomb interactions. This imposes the first require- 
ment: the singularity of the potential must be of the 
l/r-type when the distance between two particles goes 
to zero. The second requirement to the model naturally 
follows from the fact that although one cannot compute 
the potential exactly, it must have at least one "free" pa- 
rameter which describes the size of the transition region 
between l/r and Xjr^ laws (this is the parameter which 
substitutes the scattering length arising in the theory in- 
cluding only local interactions). The third requirement 
is that in the limit r — > cxd the potential must acquire the 
form (^). Finally, the potential must be differentiable in 
the whole space (except the origin) . A simple form which 
satisfies all these requirements is 
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where ad is a constant. 

A natural question which appears after the regular- 
ized potential is introduced is how it affects the results 
obtained, so far. First of all we mention that (Q) with the 
potential (^) without cutoff does not allow plane wave so- 
lutions, just because the respective integrals diverge (and 
not because of instability). Introducing a cutoff does not 
change the situation since it makes the equation explic- 
itly spatially dependent and so it still has no plane wave 
solutions. 

It is easy to see that one of the main effects of reg- 
ularization is the effect on the collapse or blow-up phe- 
nomenon. Namely, we are going to show now that any 
kind of cutoff prevents collapse. First however, we point 
out a heuristic argument in favor of our statement. The 
model (^), which is not regularized, is formally invariant 
with respect to the renormalization 5* i-^ ^'/L, r i~> Lr, 
t ^ L^t (see e.g. ||T^) while any kind of regularization 
breaks this symmetry indicating the possibility of exis- 
tence of a ground state solution. Such a solution (if any) 
will be stable if the respective Hamiltonian has a lower 
bound at constant number of particles. 

So we now proceed to prove that such a bound in the 
case of positive (or zero) scattering length exists . As 
a matter of fact after introducing a regularized poten- 
tial the proof of the stability is reduced to one given in 



pT[ . Indeed considering Eq. (|i|) and using appropriate 
dimensionless variables one can write down an associated 
Hamiltonian in the form 
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Any regularization satisfying the conditions imposed 
above satisfies also the constraint 
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where /o =const. For instance this condition is strictly 
satisfied by regularization (|l^). We would like to stress 
the essential fact that the particular function /(r) used 
to regularize is not essential provided Eq. (19) is sat- 
isfied (for example, the cut-off model also satisfies this 
requirement, as well). Next we will use the fact that 
/ y(r)|V'|2(ir > and the inequality 
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Finally 
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and thus the hamiltonian is bounded below, i.e. strict 
collapse is not possible. Of course some tendency of the 
system to compress would be observed in the dynamics. 

To conclude, in this report we have shown that the 
classical dipole-dipole interaction given by Eq. (^) is not 
consistent and must be regularized somehow to take into 
account the divergence near the origin. Fortunately the 
model is well behaved in the sense that it is possible to 
regularize in several reasonable (and physically meaning- 
ful) ways. In fact the usual numerical treatment of the 
problem includes implicitly a regularization which is why 
previous works with this interactions did not show up the 
bad-possedness of the nonregularized model. 

Finally we have shown that any reasonable choice of 
the regularization leads to a supression of the collapse 
in the sense that the ground state of Eq. ([^) exists. Of 
course the effect of other physical terms should be taken 
into consideration and depending on the scale at which 
collapse is stopped the tendency to srink the solutions 
could lead to an effective depletion of the condensate even 
in this case in which strict collapse is not possible. 
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This results means that the instability and collapse re- 
ported in and Q are related to the anomalously strong 
singularity (of 1/r^ type) of the potential at the origin. 
As we have commented before this behavior is nonphys- 
ical for small r and collapse would probably be stopped. 
This fact does not rule out however the possible existence 
of some shrinking tendency of the atomic cloud for some 
spatial scales. 
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